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ABSTRACT 

We present a two-parameter model of the extragalactic 7-ray background (EGB) in the 
0.1-100 GeV range as measured by the Large Area Telescope (LAT) onboard the Fermi 
satellite. The EGB can be fully explained as the sum of three distinct components, 
namely blazars, non-beamed AGNs (Seyfert galaxies and QSOs), and cosmic rays 
from star-forming galaxies. The contribution to the background from beamed sources 
is obtained by fitting the Fermi-hAT blazar differential number counts assuming that 
the 7-ray luminosity function is directly proportional to the radio luminosity function 
of FRI and FRII galaxies. The high energy emission from non-beamed AGNs is instead 
determined by popular synthesis models of the observed X-ray background. Finally, 
the EGB is fit by adding a third component arising from pion decay in cosmic rays, 
assuming that such component is closely linked to the cosmic star formation history. 
We find that blazars dominate at energies > 10 GeV, for E < 0.2 GeV the main 
contribution is from non-beamed AGNs, while cosmic rays are required in between. 
Because of absorption due to interaction of 7-rays with the extra-galactic background 
light, our model falls short at the highest energies probed by LAT, (> 70 GeV), 
leaving room to a possible contribution from dark matter particle annihilation. As an 
example, a particle of mass ~ 0.5 TeV and cross section (av) ~ 5 x 10^'^^ cm^ 
can accomodate the data. 

Key words: gamma rays - cosmology: diffuse radiation - BL Lacertae objects: 
general - dark matter 



1 INTRODUCTION 

The extragalactic 7-ray background (hereafter EGB) rep- 
resents a fascinating challenge since his first detection by 
SAS 2 satellite above 30 MeV (Fichtel, Simpson, & Thom- 
son 1978). The Large Area Telescope (LAT) on board the 
Fermi Gamma Ray Space Telescope has provided new EGB 
data in the range 0.1-100 GeV, after one year of observa- 
tions (Abdo et al. 2010a). The EGB is the sum of resolved 
extragalactic sources, and diffuse emission due to unresolved 
sources and/or intrinsically diffuse radiation, and its mea- 
sure requires an accurate modeling of galactic foregrounds. 
In the following we will refer to the EGB as the sum of the 
diffuse component and resolved (by FERMI) sources. The 
EGB intensity above 100 MeV as measured by Fermi-LAT 
is Jbgb(> 100 MeV) = 1.42 x lO"'^ cm"^ s"^ sr"\ where 
resolved sources accounts for ~ 27% of the emission, the rest 
being ascribed to the diffuse component. 

BL Lac objects and Flat Spectrum Radio Quasars (FS- 
RQs), collectively known as blazars, are the most common 
sources in the extragalactic 7-ray sky in the 0.1-100 GeV 
range. Therefore many authors considered these sources as 
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the main contributors to the EGB (Stecker & Salamon 
1996; Miicke & Pohl 2000; Inoue & Totani 2009; Abaza- 
jian, Blanchet, & Harding 2010a; Neronov & Semikoz 2011; 
for a review see Dermer 2007). Abdo et al. (2010b) recently 
showed that the contribution of blazars to the EGB is not 
dominant, and the focus has been shifted on different, po- 
tentially new, sources. As an example. Fields, Pavlidou, & 
Prodanovic (2010) presented a model of the cosmic produc- 
tion of 7-rays powered by the interaction of cosmic rays with 
the interstellar medium in star-forming galaxies. The model 
reproduces quite well the behavior of the EGB up to 10 GeV, 
but it is not fully satisfactory at higher energies. 

In the present letter, we propose a model of the EGB 
based on the 7-ray emission from blazars, non-beamed 
AGNs (Seyfert galaxies and QSOs), and star-forming galax- 
ies. We show that a combination of these three populations 
fully explains the EGB in the Fermi-LAT energy range. 



2 EXTRAGALACTIC 7-RAY SOURCES 

Our model employs three components, blazars, AGNs, and 
star-forming galaxies. First, we assume that the blazar 7- 
ray luminosity function (LF) is proportional to the radio LF 
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of FRI and FRII galaxies. To convert radio to 7-ray lumi- 
nosities we use the spectral energy distribution proposed by 
Fossati et al. (1998). AGNs are taken from population syn- 
thesis models of the X-ray background. For the star-forming 
galaxy emission we adopt a model of pion decay in cosmic 
rays. 

Our strategy is the following: we first fix the contribu- 
tion of blaaaxs to the EGB by fitting the Fermi-LAT differ- 
ential number counts, having as sole fitting parameter the 

fraction of radio galaxies beamed toward the Earth, hence 
appearing as blazars. Then we add the fixed background 
from non-beamed AGNs (Inoue, Totani, & Ueda 2008). Fi- 
nallj', the contribution of galaxies is obtained by fitting the 
EGB. Note that since the shape of the cosmic ray contribu- 
tion is determined by pion decay and star formation history, 
the only free parameter is the star formation "efficiency" of 
molecular hydrogen (Stecker & Venters 2011). We checked 
that results do not change if we instead perform a simul- 
taneous, 2-parameter fit to the number counts and to the 
background data. 

In the following subsections we briefly describe the de- 
tails of our modeling. Throughout the paper we consider a 
concordance cosmology with Ho = 70 km/s/Mpc, Orn = 0.3, 
and = 0.7. 



2.1 Blazars 

The blazar contribution (in photons s~^ cm~^ sr~^ MeV~^) 
to the EGB at the observed energy Eq is 
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where d$-y [L^,z)/d log L^, is the 7-ray LF and is vL^ (in 
erg/s) at 100 MeV, dn{L^, z)/dE is the unabsorbed photon 
flux per unit energy E = Eo{l -\- z) measured on Earth of 
a blazar with luminosity at rcdshift z, and t^~.i[Eq,z) 
is the optical depth for 7 — 7 absorption. We adopt the 
extragalactic background light (EBL) model by Finke, Raz- 
zaque, & Dermer (2010). In the above equation dV/dz is 
the comoving cosmological volume. We set logL™'" = 43 
and log L!^"" = 50. 

The number of sources A''(> Fph) per steradian with 
photon flux greater than Fph is 
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The 7-ray LF of blazars is presently uncertain (for an esti- 
mate see, e.g., Abdo et al. 2009), so that one has to rely on 
the LFs computed in other bands, e.g.. X-rays (Narumoto & 
Totani 2006; Inoue & Totani 2009; Abazajian et al. 2010a), 
or radio (Draper & Ballantyne 2009). We adopt here the 
radio LF at 151 MHz of FRI and FRII (Willot et al. 2001), 
assumed to be the parent populations of blazars: 
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where Lr is uL^ at 151 MHz, and the constant k is the 
fraction of blazars over all radio galaxies, and it is out flrst 



fit parameter. In order to convert radio into 7-ray luminos- 
ity, wc rrmst rely on the blazar spectral energy distribution 
(SED). According to the SSC model, the blazar SED shows 
two broad peaks, the first located in the IR/X-ray band, due 
to synchrotron emission from the rclativistic jet, and the 
second in the 7-ray band, due to the inverse Compton scat- 
tering of synchrotron photons (e.g., Ghisellini et al. 1998). 
As pointed out by Fossati et al. (1998), a relation between 
the radio and/or bolometric luminosity and the energy of 
the two peaks exists (the so-called blazar sequence, sec also 
Donato et al. 2001). We use the SEDs computed by Inoue 
& Totani (2009) based on the empirical determinations of 
Donato et al. (2001). 



2.2 Non-beamed AGNs 

Hot thermal electrons in accretion disk coronae of Seyfert 
galaxies and QSOs scatter off UV photons into the hard 
X-ray/7-ray bands (e.g., Haardt & Maraschi 1991, 1993). 
Such emission, filtered by absorbing material at a distance 
of several parsccs form the central source, is though to be 
responsible of the observed X-ray background (e.g., Sctti 
& Woltjer 1989; Madau, Ghisellini, & Fabian 1994; Gilli, 
Comastri, & Hasinger 2007) . According to Inoue et al. (2008) 
the AGNs responsible of the hard X-ray background can also 
account for the 1-10 MeV EGB, provided that an extra non- 
thermal electron population is considered. As in Abazajian 
et al. (2010a), we use the model of Inoue et al. 2008 with a 
powerlaw index of non-thermal electrons = 3.5, a lower limit 
of the Lorentz factor distribution = 4.4. 



2.3 Star-forming GalcLxies 

The 7-ray spectrum of a star-forming galaxy is based on 
the assumption that 7-ray emission is due to the decay of 
tt" mesons. The tt" mesons form in the inelastic collision be- 
tween cosmic rays and the ISM. According to Stecker & Ven- 
ters (2011), the specific 7-ray photon spectrum Lph (photons 
s~^ MeV~^) of a star- forming galaxy is related to the av- 
erage pionic 7-ray production spectrum per hydrogen atom 
(qMEo)) (Dermer 1986; Mori 1997) as. 



iph(So) = (qn{Eo))Nn, 



(4) 



where A'h is the total number of hydrogen atoms in the 

galaxy, both in atomic and molecular form. 

We adopt the Strong Coupling 7-roj/ - Star Formation 
Rate Model of Stecker & Venters (2011), where A^h is related 
to the star formation rate. According to the model, the star- 
forming galaxy contribution to the 7-ray background is 



^ ,p , 1 r , di (l + 7^) , . , , _ 

where dl/dz — cH~^{l + z)~^ with H{z) the Hubble param- 
eter, TZ ~ 0.9 is the ratio of atomic-to-molecular hydrogen 
density in star- forming galaxies (see Leroy et al. 2008), and 
PSFR is the cosmic star formation rate (we use the fit pro- 
posed by Li 2008). The parameter ^(-^2) (the star forma- 
tion efficiency of molecular hydrogen, sec Bigiel et al. 2008, 
Gnedin et al. 2009; Bauremeister et al. 2010) is the ratio be- 
tween psFR and the cosmic density of molecular hydrogen, 
and it is our second fit parameter. 
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Figure 1. The data points are the differential blazar number 
counts detected by Fermi-hAT as a function of the integrated 
flux above 100 MeV Fioo (Abdo et al 2010b). The solid red line 
is our best fit. 
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Figure 2. Upper panel: the EGB measured by Fermi-hAT (data 
points by Abdo et al. 2010a), and our fit to it. The total emission 
(solid black line) is shown together with the contribution from 
blazars {long-short dashed blue line), non-beamed AGNs {short 
dashed red line), and 7-ray emission from cosmic rays in star- 
forming galaxies {long dashed green line). Lower panel: fractional 
contributions of the three different components to the EGB. 



3 RESULTS 

We use the model discussed in the previous section to fit 
the Fermi-LAT EGB. As already discussed, we first com- 
pute the contribution of blazars by fitting their Fermi-ljAT 
differential number counts (Abdo et al. 2010b), then we add 
the background from non-beamed AGNs (Inoue et al. 2008), 
and finally we compute the contribution of galaxies by fit- 
ting the EGB. We obtain the same best fit values (within 
the errors) if we perform a simultaneous fit to the number 
counts of blazars and to the background. Results and num- 
bers reported here refer to this second approach. We remark 
that the shape of the star-forming galaxies and blazar com- 
ponents are determined a priori, while only their normaliza- 
tions are allowed to vary. The contribution of non-beamed 
AGNs is instead totally fixed. 

In fig. [T] we show our best fit to the differential blazar 
number counts, while the best fit to the 7-ray background 
is shown in fig. (2] upper panel. In the lower panel of fig. [2] 
we show the relative fractional contributions of the three 
different components. 

The best fit parameter values we obtain are k = 
3.93 ± 0.01 X 10"" and ^(^2) = 4.07 ± 0.4 x 10"^° yr~\ 
with xVd-o.f = 1.15 for 29 d.o.f.. The best fit value of 
^, which represents the star formation efficiency of molec- 
ular hydrogen is in agreement with observational values 
(C ~ (5.25 ±2.5) X 10""'yr"\ Leroy et al. 2008). The num- 
ber ratio of blazars to radio galaxies k can be thought as a 
measure of the beaming factor of the relativistic jet, which in 
turn is related to the bulk Lorentz factor F. From k ~ 1/2F^ 
we derive F ~ 35. 



4 DISCUSSION AND CONCLUSIONS 

We proposed a two parameter model for the extragalactic 7- 
ray background (EGB) measured by Fermi-LAT. We showed 
that the EGB can be explained by a combination of blazars, 
non-beamed AGNs (Seyfert galaxies and QSOs), and high 
energy photons originating by the interaction of cosmic rays 
and ISM in star-forming galaxies. 

Our EGB model differs with respect to other existing 
similar studies (e.g., Inoue & Totani 2009; Abazajian et al. 
2010a; Stecker & Venters 2011) in the 7-ray LF of blazars 
that we derive directly from the radio LF allowing only a free 
relative normalization (i.e., the LF faint and bright ends are 
fixed before fitting), and/or in the use of the Fossati et al. 
(1998) SED instead of a simple power law parametrization 
of blazar high energy spectrum. 

We found that non-beamed AGNs dominate the EGB 
below 0.2 GeV, while blazars take over above 10 GeV. Cos- 
mic ray emission is then needed at intermediate energies at a 
level ~ 50% of the total. Our model parameters are the rela- 
tive ratio of blazars to radio galaxies, and the absolute nor- 
malization of the contribution from star-forming galaxies. 
The best fit value of the relative number of blazars with re- 
spect to radio galaxies can be translated into a bulk Lorentz 
factor of the relativistic jet F ~ 35, larger than the average 
value F ~ 15 estimated by Ghisellini et al. (2010). The two 
values could be reconciled if blazars commonly show secular 
7-ray large variability which modulates the 1-year average 
fiux, as recently proposed by Ghirlanda et al. 2011. Our 
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Figure 3. Upper limits of the cross section {av) as a function of 
particle mass rriDM for annihilating DM. The lower (upper) curve 
represents the l(2)-tT limit. See text for details. 



model fit to the EGB constrains the so-called " star for- 
mation efficiency of molecular hydrogen", which we found 
well within existing, much looser observational constraints 
(Leroy et al. 2008). 

Clearly, the specific best fit values obtained depend 
upon the details of our model, in terms of star formation 
rate adopted, models for the 7-ray emission of star-forming 
galaxies, blazar LF and SED. Nevertheless the overall pic- 
ture appears quite robust, with an important role played by 
star-forming galaxies, with blazars dominating only at the 
higher energies probed by Fermi-LAT. 

Though our model is statistically fully acceptable, it is 
interesting to note that the highest data point of the EGB 
(see fig. 2 , upper panel) lies above our best fit model. In 
the energy band 50-100 GeV absorption of 7-rays due to 
the interaction with the EBL is significant. Different theo- 
retical EBL models have been proposed in the last few years 
(e.g., Franceschini, Rodighiero, & Vaccari 2008; Gilmore et 
al. 2009; Kneiske & Dole 2010; Finke et al. 2010; Dommguez 
et al. 2011), resulting in somewhat different optical depth 
for photon-photon interaction. As already discussed, we fol- 
low Finke et al. (2010), and we checked that even adopting 
the model of Kneiske & Dole (2010), which gives the low- 
est 7-ray absorption, our EGB model still falls short in the 
70-100 GeV range. A possible, intriguing explanation is the 
presence of an extra emission from annihilating dark mat- 
ter (DM) particles (see, e.g., UUio et al. 2002). Recently, 
Abazajian, Blanchet, & Harding (2010b) performed a de- 
tailed analysis of possible DM candidates in the context of 
Fermi-LAT EGB. A full discussion of DM particle contribu- 
tion to the EGB is beyond the scope of the present Letter. 
For illustrative discussion, here we adopt a specific annihi- 
lating DM model, and compute its contribution to the EGB. 

The 7-ray background produced by annihilating DM 
is calculated following Abazajian et al. (2010b) and Ando 



(2005). We found, as an example, that a particle of mass 
~ 0.5 TeV and cross section {av) ~ 5 x 10~^^ cm"^ s"""^ can 
easly accomodate the last data point. However its presence 
is not statistically required by the fit, so it is fair to consider 
only upper limits to the DM component. Fig.[3]shows our re- 
sults in terms of cross section (av) and particle mass moM- 
The lower (upper) curve is computed by adding the DM 
background to our EGB model, allowing a increase of 1 
(4) with respect to the best fit, hence representing the l(2)-(j 
upper limits of {av) for a given muM- As an example, as- 
suming {av) = 3 X 10^^® as required for leaving the observed 
relic density of DM (Jungman et al. 1996), we can exclude at 
l(2)-cr level DM particles with moM < 100(10)GeV. More 
massive particles can have a larger cross section, and still be 
compatible with EGB data. Our limits are consistent with 
other, more refined, determinations (e.g., Abazajian et al. 
2010b). 

Finally, it is worth asking whether DM particle annihi- 
lation signal could substitute the 7-ray emission from star- 
forming galaxies in fitting the EGB. We checked that the 
joint fit to blazar number counts and EGB adopting an al- 
ternative three component model (i.e., blazars, non-beamed 
AGNs, and DM) is statistically unacceptable (x^ ~ 47.5 for 
28 d.o.f., to be compared to = 33.3 for 29 d.o.f. as dis- 
cussed in section 3). We conclude that star- forming galaxies 
seem to be a necessary, and important, component of the 
EGB, while DM is not strictly required. 
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